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X/^ \ A crucial element of scattering theory and the LSZ reduction formula is the 

00 ' assumption that the coupling vanishes at large times. This is known not to 

hold for the theories of the Standard Model and in general such asymptotic 
dynamics is not well understood. We give a description of asymptotic dynamics in 
field theories which incorporates the important features of weak convergence and 
physical boundary conditions. Applications to theories with three and four point 
interactions are presented and the results are shown to be completely consistent 
' with the results of perturbation theory. 
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X 1 Introduction 

Descriptions of scattering in quantum field theory assume that at large times the par- 
ticles are widely separated and behave like free particles. This assumption, which 
underlies the LSZ formalism, is incorrect for many theories. The most obvious exam- 
ple of this is when the incoming or outgoing system includes bound states, but it also 
fails if the physics is characterised by long range interactions. Since most of the physics 
of the standard model falls into at least one of these categories (confined quarks, mass- 
less gauge bosons) it is very important to have a precise understanding of the dynamics 
of quantum field theories at large times. 

Generally, then, it is assumed that at asymptotic times the Heisenberg fields become 
free ones 

lim (j)(x) - ZV^outCx) , (1) 

t— >oo 
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and similarly for t — > — oo. We should note that this behaviour can only be taken to 
hold as a weak limit, between matrix elements, since otherwise (see Sect. 5-1-2 of [|I]]) 
one can show, from the Kallen-Lehmann representation, that the fields are free at all 
times. 

The limit in ([!]) is sometimes discussed in the framework of an 'adiabatic approxi- 
mation', in which the coupling constant is taken to be multiplied by a function which 
is one during the scattering process and approaches zero for very large (positive or neg- 
ative) times. This is unsatisfactory since it assumes the desired answer which ought 
rather to emerge from the theory itself. It can also be wrong, as in the case of Quantum 
Electrodynamics (QED). 

QED, the paradigm for the Standard Model, has long range interactions. The 
masslessness of the photon means that the potential between static charges falls off 
only as 1/r. It has been known for a long time p],|3| that this means that ([!]) does 
not hold and that any attempt to impose such a relation generates infra-red (IR) 
divergences in the wave-function renormalisation constant of (charged) matter fields. 

This has been studied M in the relativistic theory by Kulish and Faddeev (KF) and 
their general approach to asymptotic dynamics has been utilised by various authors, 
see, e.g., [PHT3|- We shall now give a brief sketch of the procedure adopted by KF and 
what their results seem to indicate for QED. 

They considered the usual QED interaction Hamiltonian 



Hmt (t) = -eJd 3 xA^t,x)J^t,x), (2) 

where J M (t, x) = ip(t,x)'y ll 'if)(t,x) is the (conserved) matter current. In order to carry 
out the LSZ reduction of the S'-matrix we must be in the interaction picture so that, 
although the time evolution of the states is determined by @, the evolution of the 
fields themselves is given by the free Hamiltonian. One may then insert the free field 
expansions in (H). These plane wave expansions are 

= Jj0js ^= {HP, s)u s { P )e-^ + dt(p, s)v s (p)e^} (3) 

where the notation implies a sum over the s indices. Working in Feynman gauge we 
have 

Mx) = Jwfir k i a ^ e ~ lkx + a U k y k ' x } > ( 4 ) 

where ho = E p = a/|j>| 2 + rn? and uJk = \k\. Inserting these into (0) results in eight 
terms which may be grouped according to the positive and negative frequency compo- 
nents of the fields. Each of these pieces will have a time dependence of the form e % ^ 1 
where ip involves sums and differences of energy terms. 

KF claimed that, for \t\ — > oo, only terms with ip tending to zero contribute to the 
asymptotic dynamics. Since the spatial integration in (|2]) generates a momentum delta 
function, only four terms with ip = ±(E p+ k — E p ± u>k) would then have a large t-limit. 
This vanishing of E p+ k — E p ±Uk ~ can only take place in QED because the photon 
is massless, and it only occurs for soft photons, i.e., for Uk ~ 0. This is in accord 
with perturbation theory: the breakdown of the S'-matrix occurs for soft photons and 
giving the photon a small mass acts as a cut-off on these divergences. An asymptotic 
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approximation to (|2|) is obtained from the lowest order term of the Taylor expansion, 
in powers of k, of the Hamiltonian. This yields 




e / d 3 x A^(t,x)J^(t,x) 



(5) 



where 




d 3 p p^ 



p(p)S 3 (x 



P_ 



v 




(6) 



and p(p) is the charge density 



P{p) = Yl s ) & ^' s ) ~ d ^P> s ) rf ^' s ) 



(7) 
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According to KF therefore, the asymptotic Hamiltonian is the integral over all momenta 
of the current associated with a charged particle of velocity p^ / E p . This non- vanishing 
Hamiltonian finds its perturbative expression in the branch cuts, rather than poles, in 
the on-shell Green's functions in the matter fields of QED. An attractive aspect of this 
theory is that it completely dispenses with the adiabatic approximation. Although this 
discussion seems to pick up the problems in applying the LSZ scheme to QED and, in 
particular, correctly identifies the problem with long wavelength (ojk ~ 0) photons, it 
cannot be regarded as the end of the story. 

These arguments are used extensively in other theories, such as QCD, where the 
physics is not well understood, and where greater reliance is put on the mathematics. 
On the other hand the theory is employed at the level of operators whereas it is more 
appropriate, in quantum field theories, to work at the level of matrix elements and weak 
limits. The KF approach also makes no connection between the large time limit and the 
separation of particles at large distance. As we shall see below, the naive application 
of this approach to massive 4 theory would indicate that the LSZ formalism should 
also break down, which it very evidently does not. 

This paper is concerned with constructing a new approach to asymptotic dynamics 
within the context of weak convergence, and with appropriate physical boundary con- 
ditions corresponding to the separation of particles. We will apply it to a variety of 
interactions and show that it yields results which are completely consistent with what 
is known from explicit perturbative calculations. 

In Sect. 2 we will show that the KF argument cannot be applied to four point inter- 
actions. Sect. 3, which is the heart of this paper, develops the method for massive 4 
theory and then applies it to both the three and four point interactions of scalar QED. 
In this way we will see both the well known spin independence of the IR structure 
in the abelian theory with massive charges and, for the three point vertex, regain the 
results of |§]. A discussion of the implications of these results for perturbative calcu- 
lations in the standard model is presented in Sect. 4. Some technical details are given 
in appendices. 

2 Four Point Interactions 

As we have indicated in the introduction, there are objections to the KF view of 
asymptotic dynamics: firstly, that the statements have been framed specifically for the 
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picture of the u k ~ infra-red theory, which is well known to suffer from divergences, 
and not as a general statement on the behaviour of asymptotic limits; secondly, and 
more critically, that the prescription does not translate to other quantum field theories. 
It is the latter objection which is the most serious and the one which we shall now 
demonstrate. 

To be precise, we shall show that the KF argument applied to the case of massive 4 
theory is not sufficient to show that the asymptotic limit of the coupling term vanishes. 
However, 4 theory is the standard textbook example of an interacting quantum field 
theory and its perturbation theory is straightforward - the coupling must vanish for 
well separated particles . 

To begin then, let us take the standard free field expansion of the scalar field 

d 3 k 



with E k = a/I^I 2 + m 2 and m the mass of the particle. The commutator relations are 

[a(k) , a\k')] = (2tt) 3 2E k 5 3 {k - k') (9) 
and the interaction part of the Hamiltonian for the theory under consideration is then 

Hmt = ^. J d3x : ^ (a;,t) : ' (10) 

where the : : indicates normal ordering. 

When the expansion (H) is inserted into (|T0|) and the resulting expression is sim- 
plified, then, after normal ordering, it will be found to consist of twelve terms, each 
of which has an exponential term where the exponent is made up of sums and differ- 
ences of the energy eigenvalues. Some of these exponents are obviously non-vanishing. 
According to the methods of KF which were described in the introduction, the inte- 
grals containing these exponentials may be ignored. However, not all of the integrals 
involved have exponents which are so easily dealt with and one of these, which we shall 
now consider, is 

r d 3 pd 3 qd 3 k ot(fe) a t(p + q - k)a(p)a(q) it(EB+Ea _ Ek _ En+a _ k) 
J (2tt) 9 2E p 2E q 2E k 2E p+q _ k ' ■ K ) 

Applying the methods described in the introduction, any non-vanishing asymptotic 
dynamics will come from those terms and those momenta for which the exponent 
vanishes. To this end we must determine if the equation 

E p + E q — E k — E p+q _ k = (12) 

has any solutions. 

Far from this being a difficult problem, upon reflection it becomes obvious that there 
are infinitely many solutions! This is most easily seen by noting that the problem is 
equivalent to that of finding solutions to the simultaneous system of equations 



\/|p| 2 + m 2 + a/I^I 2 + m2 = \/|fc| 2 + m2 + V\l\ 2 + m2 

p + q = k + l. (13) 
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Figure 1: Triangle of Momenta 



The incoming and outgoing momenta must have some connection with a scattering 
process so that the obvious, trivial solutions, which can be found by taking, say, q = 
1 = 0, will be ignored. 

Referring to Fig. [I], let p =AB, and q =BC . Now imagine pivoting the rigid 
triangle ABC, with the line AC as the hinge, to get a new triangle which is congruent 

to ABC . If the new triangle is AB'C then take k =AB' and / =B'C . The vectors 
p, q, k, I then automatically satisfy both the conditions of ([13]). 

One can do this with any two non-parallel vectors, for example (with m = 1) take 
p = (l,l,0),q = (l,-l,0),fc = (l,l/V2,l/V2),l = (1, -1/V2, -1/V2) • Another 
example is p= (-4/v^3, 2V2/V5, 0), q = 0,k = (-2/v^, v^/v^, 1), and I = p+q-k . 

No physical meaning could be given to such an arbitrary set of momenta. Following 
the arguments of KF however, we conclude that the existence of this wide range of 
momenta for which (12) vanishes, indicates that the asymptotic dynamics of the system 
is determined by this set of points. This would suggest that problems could arise in 
the associated perturbation theory for massive 4 theory, and this is well known to be 
false for this textbook example of a quantum field theory. 

At this point we wish to stress that this apparent contradiction has implications 
beyond this simple example as there are many important quantum field theories, such 
as QCD and the Higgs mechanism, in which four point interactions have a major role. 
If there is a deficiency in the understanding of the asymptotic dynamics in this simple 
scalar four point theory then it is difficult to see how one might proceed on this basis, 
with any confidence, in other theories of the standard model. It is clear from this that 
the KF argument needs considerable refinement if it is to be applied to the standard 
model. 



3 General Approach to Asymptotic Dynamics 

The major flaw in the KF argument is its reliance on strong operator convergence. 
Instead of concentrating on the Hamiltonian itself, we shall focus our attention on 
matrix elements. The limits we shall consider are weak limits and this is in keeping 
with the LSZ formalism in quantum field theories. 

If ^in represents an incoming wave packet in this scattering experiment and \I/out is 
an outgoing wave packet, then the matrix element of interest is <\E f Q UT |7i int |\I / IN >. This 
is a time dependent c-number and the elementary notion of convergence may be adopted 
in the investigation of its asymptotic limit. In the case of massive scalar 4 it will be 
shown that, under conditions which have a straightforward physical interpretation, its 
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asymptotic limit is zero. 

This section is organised as follows. In Sect. |3.1| we shall use massive 4 theory 



as a vehicle to study asymptotic dynamics. Having established the basic method and 
shown that it is consistent with what is already known for 4 theory, we shall then 
employ it to determine the asymptotic dynamics of scalar QED. 

There is a particularly interesting feature of scalar QED which makes the study of its 
asymptotic dynamics worthwhile. The interaction Hamiltonian in scalar QED consists 
of a sum of two terms, a three point interaction term, similar to that in (fermionic) 
QED, and a four point interaction term, which has no parallel in QED, so that the 
dynamics of scalar QED is richer than that of fermionic QED. From perturbation 
theory it is known that (fermionic) QED and scalar QED have the same infra-red 
problems, so their asymptotic dynamics are the same. Perturbation theory also tells 
us that the four point interaction term is divergence free so it must have vanishing 
asymptotic dynamics. We would expect this to emerge, in a natural manner, from any 
satisfactory theory which describes the asymptotic dynamics of scalar QED. Further, 
since the infra-red problem is spin independent, the asymptotic dynamics of the three 
point interaction will, with the obvious changes, be the same in either fermionic or 
scalar QED. 

We shall begin our examination of scalar QED in Sect. |3.2| with the quartic term, 
which turns out to be the easier to deal with. We shall show that its asymptotic 
dynamics is similar to that of the scalar 4 theory and has a zero limit. After this, in 
Sect. |3l| we shall turn our attention to the cubic term and the infra-red problem. We 
shall prove that the asymptotic dynamics for the infra-red problem is exactly the same 
as that for a system in which the Hamiltonian is derived from a current associated with 
a moving charged particle with known, non-trivial asymptotic dynamics. 

3.1 Scalar 4 Theory 

Consider the following incoming and outgoing wave packets. 

^i N = J d 3 rd 3 w f(r)g(w)a ji (r)a^{w)\0> , 



\1>out = J d 3 ud 3 v h(u)i(v)a)(u)a){v)\{]> 



(14) 



with the functions /, g, h, i being test functions for the wave packets. Referring to (§), 
@, (|T0D , and the expressions (0), one finds that <^o UT |7ii nt |\I/iN> reduces to a single 
integral which is proportional to 

J d 3 pd 3 qd 3 kh{k)i{p + q - k)f{p)g{q)e- iti) , (15) 

with the exponent ip having the value i/j = E p + E q — — E p+q _k . 

Notice that the exponent in this term has essentially the same structure as the 
term in ( fTT| ) which appeared to cause problems when applying the methods of KF. 
The difference now is that ( |15|) is a straightforward integral and not an operator, so 
that elementary methods may be applied to find the asymptotic limit. The machinery 
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we shall employ is the method of stationary phase JT4|. Briefly, this says that, provided 



there is no point in the region of integration at which all of the first order partial 
derivatives of ip are zero, then the integral ([15]) tends to zero as \t\ — > oo . 

The terms in ijj have the form E\ = ^/\l\ 2 + m? and, since m/0, they will all have 
first-order partial derivatives for all values of I. The first order derivatives of ip are then 
given by 



dtp 


Pi_ 


Pi + qi- h 


dpi 


E p 


Ep+q—k 




1j 


Pj + Qj ~ k j 


d qj 


E q 


Ep-^q—k 


d^jj 


k n 


Pn ~\~ Qn k 


dk n 


Ek 


Ep+q—k 



(16) 



If at some point all of these are zero then, in particular, dip/dpi = dtp /dqj = 
for all possible values of i,j . This implies that p/E p = (p + q — k)/E p+ q_ k and 
q/ E q = (p + q — k)/E p+q _ k so we must also have p/E p = q/E q . If we can exclude the 
set of points for which this condition holds from the domain of integration, then the 
integral in (|TJ) will vanish as \t\ — ► oo . 

In (|15|) , the test functions /, g for the the incoming wave packet have the variables 
p, q as arguments. The expressions p/E p and q/E q represent the velocities of the 
respective incoming wave packet. Experimentally, scattering is prepared by setting up 
the apparatus in such a way that the two beams of particles are brought together from 
different directions, i.e., with different velocities. This information can be incorporated 
into the incoming wave packet by ensuring that the supports of the test functions 
exclude the possibility that p/E p = q/E q . The precise statement of the requirement is 
that the test functions f,g must have non- overlapping supports in velocity space. This 
condition on the test functions, of having non-overlapping supports in velocity space, 



is central to the construction of the S'-matrix (see Sect. 13.4 of [[TBI ). 

To restate, if the test functions /, g have non-overlapping supports in velocity space 
then the integral in ( |i~5"|) vanishes as \t\ — > oo . This is exactly the behaviour that 
one would expect for this particular scattering process but a further question remains 
to be answered: what constraints does this choice of test functions for the incoming 
wave packet impose on the outgoing wave packet? If this picture is to display all of 
the features of this particular scattering process then some combinations of outgoing 
particles must be excluded. The outgoing particles must behave as free particles at 
asymptotically large times which means that, as was the case for the incoming wave 
packets, their test functions must also have non-overlapping support in velocity space. 
However, this must be a consequence of the condition imposed on the test functions 
for the incoming wave packet and not an independently imposed condition. 

The arguments of the test functions for the outgoing wave packet are p + q — k and 
k, and the equality of these two variables is equivalent to p+q = 2k. This is simply the 
expression of conservation of momentum. Another principle in any scattering theory is 
conservation of energy. In this case this is expressed as E p + E q = 2E k . Finally then, 
we must show that the two conditions 

p + q = 2k, 
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E p + E q = 2E k , (17) 

are incompatible with the functions /, g having non-overlapping support in velocity 
space. Since the masses of the two incoming particles are equal, the equation p/E p = 
q/E q is equivalent to p = q and non-overlapping in velocity space is equivalent to 
non-overlapping in momentum space. In this case we need to show that if p ^ q then 
the conditions of flTTD are impossible. 

Again, for simplicity, let us take the mass m = 1 . The first equation in fll7f) means 
that the vectors p, q, k, are coplanar. In that case, we may choose a unit vector n 
orthogonal to this plane and write p' = p + n, q' = q + n, k' = k + n . Since n is 
orthogonal to the plane of p, we have \p'\ = a/|p| 2 + 1 = E p , with similar expressions 
for q, k . This means that fllTf ) may now be written as 

p' + q' = 2k' 
\p'\ + \q'\ = 2\k'\. (18) 

From the triangle inequality we know that this is only possible when p' and q' are 
parallel, i.e. if there is a number A such that p' = Xq' . In terms of p, q and n, this 
can be rearranged into the form p — Xq = (1 — A)n . The only solution for this is with 
A = 1 so that p = q . 

There are other possible matrix elements that are associated with the four point 
interaction term. There is the possibility that the incoming wave packet consists of a 
single field, with the outgoing wave packet made up of three fields, and there is the 
reverse. Both of these cases may be treated in precisely the same manner as the above 
case, and with precisely the same conclusions. We omit the details. 



3.2 Scalar QED: the four point interaction. 

The interaction in scalar QED is more complicated than that of QED due to the 
existence of the extra term representing a four point interaction. In this section we 
shall study the asymptotic properties associated with this term and show that it has 
trivial asymptotic dynamics. This will require extending our techniques since we shall 
have to deal with wave packets which have massless particles as an essential part of 
their structure. 

The method of stationary phase that was used to determine the limit of flT5|) , in 



Sect. |3TI|, was dependent upon some of the properties of the partial derivatives ([LED °f 
the various energy eigenvalues. In the case of massless particles the energy eigenvalues, 
which are of the form Uk = \k\, are not differentiable at the origin but this will not be 
a barrier to the application of this technique. 

We begin by writing out the full, normal ordered interaction Hamiltonian for scalar 
QED, which is 

n int (t) = -e [ d 3 x : J»{x)AJx) : (19) 



with : : being normal ordering and where the current J M is given by 

= U?-eJ% (20) 
= i ( + Jf 2 ) — e J% , 
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with the obvious meaning given to the components of J M defined in (|20|). 

We shall work in Feynman gauge and take the plane wave expansions given by 



= 0+0) + <f>-(x) 



(21) 



/d 6 k 1 

= A+(«)+i4;(a!). 
The commutator for the photon is given by 

K(p), al(q)] = -{2iif2u k g^\p - q) . (22) 
When (plD is substituted into (|19"D and rearranged, the quartic interaction will be 



found to consist of 12 terms, each of them with differing exponential terms corre- 
sponding to the different possible interactions. We shall consider the matrix element 
represented by the following diagram. 




We shall take as our wave packets the following expressions 

#in = Jd 3 rd s w f{r)tf(r)(?{w)al{w)\0>, 

(23) 

#out = Jd 3 ud 3 v g(u)tf(u)h»(v)al(v)\0> , 

where the /, cf, g, are the respective test functions. One then finds that the ampli- 
tude <^out I ?^int(£)| ^in> is a single term, constructed from J% , which is an integral 
proportional to 

J d 3 pd 3 qd i kc\p)h^q)g{k)}\k + q-p)e^\ (24) 

where ijj = u p + E k+q _ p - E k -u q . 

Notice that the Einstein summation convention means that this integral is actually 
a sum of three integrals, with i = 1,2,3. The incoming and outgoing charged fields 
must be separated, i.e., the test functions / and g must have disjoint support. If not, 
then, by conservation of momentum, there will be no separation of the incoming and 
outgoing photons. In that case, no scattering will have taken place. Thus it must follow 
that, for each i, the test functions q, hi must have disjoint support so that at least one 
of them will not have the zero vector in its support. Without loss of generality, let us 
suppose that this is q. 
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The function ip will then have continuous partial derivatives in p, k and these are 
given by 



dip Vi_ _ kj + Qi- Pi 

dpi UJ p Ek+ q -p 

(25) 

dip kj + qj — pj kj 
dkj Ek+q-p Ek 

The vanishing of these expressions for all i, j would then imply that p/w p = k/ Ek which 
is impossible, since the former is a unit vector while the latter is not. The method of 
stationary phase (applied only to the variables p, k) can again be used to prove that 



(|24| ) will vanish as \t\ — > oo . 

There are other possible scattering events that are described by this four point 
interaction Hamiltonian, e.g., when the wave packets consist of two incoming photons 
and two outgoing charges (or vice versa) . If the photons are separated according to our 
scheme (their test functions have disjoint support), then a similar exercise will show 
that this picture also gives rise to vanishing asymptotic dynamics. 

3.3 Scalar QED: the infra-red approximation 

The cubic term in the interaction Hamiltonian (|19D, which comes from the term in 
fl20D , is made up of two parts, J{\ and Jf 2 . These two terms are similar in their structure 
and both contribute to the problem of infra-red divergences. We shall consider only 
the first of these, Jf 1; the results for Jf 2 being substantially the same. 

As in QED, the infra-red problem in scalar QED occurs in relation to a scattering 
process in which an incoming charged particle emits a photon. 




We shall consider the case when the wave packets are given by 

*in = Jd 3 yf(y)tf(y)\0> 

^out = Jd 3 ud 3 v g(up(u)h fM (v)al(v)\0> . (26) 

The matrix element will then be found to be given by 

<*J> UT | H int (t)\ * IN >= -e J d 3 qd 3 kf(q + k)g(q)q^h^k)e-^ (27) 

where ip = E q+ k — E q — Uk . Now it is easy to see why the previous method cannot be 
applied in this case. Since the photon is massless, the corresponding energy eigenvalue 
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is Wfe = |fe|. The expression for if), therefore, will not have partial derivatives in fc, at 
\k\ = 0, and the partial derivatives in qj will vanish when \k\ = 0, for any value of q. 
It is this regime, when \k\ =0, which gives rise to the problem of infra-red divergences 
and the difficulties associated with its asymptotic dynamics cannot be avoided. 

Following KF we shall compare the asymptotic dynamics of this scattering process 
with that governed by the scalar version of the asymptotic Hamiltonian (H). 

The system that we shall take for our comparison is then the one defined the 
asymptotic current given by 

JM = 1 1 wf Gk) 2 6,(9)6(9) ^ 3( * - ^ • < 28 » 

where we omit terms which do not contribute to this matrix element. The interaction 
Hamiltonian is given by 

WS = -c J d 3 x J^ s (x) A„(x) . (29) 

If we take the incoming and outgoing wave packets (^6|) and contract, in the usual 
fashion, with 7i^ t (t) then we obtain the amplitude 

<*&ur I *in>= -e J d 3 q d 3 kf(q)g(q) q ^h fl (k)e-^' t } (30) 

where ip' — q ■ k/E q — Uk ■ 

The asymptotic dynamics of the systems defined by the interaction Hamiltonians 



(ri9| ) and (|29|), and the wave packets (0), will then be the same if we can prove that the 
difference between the integrals in ([27]) and (|30D vanish for asymptotically large time. 
This amounts to showing the following theorem whose proof is given in the Appendix. 



Theorem 1 Let f,g,h^,ip,ip' be defined as in ( \2Tj ) and (pfy). Then 



lim / d 3 qd 3 k[f(q + k) e ^ t -f(q) e ~^ t ]g(q)q^h fl (k) = 0. (31) 

This shows that the physics of scattering in scalar QED at large times is fully described 
by the asymptotic Hamiltonian ( p9|) . The four point interaction vanishes and the three 
point one can be described using the simple current (p8|). We have thus extended 
the results of KF to scalar QED and proved the spin independence of the asymptotic 
dynamics in abelian gauge theories. 



4 Conclusions 

It is not necessary to assume that the coupling constant asymptotically switches off. 
As we have seen, one can, for theories like massive 4 , prove that the asymptotic dy- 
namics is free or, for theories like QED, with rather more effort, determine the form of 
this asymptotic interaction. The arguments for determining the asymptotic properties 
of interactions in quantum field theories, proposed by Kulish and Fadeev have been 
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improved upon and made applicable to a more general type of interaction, including 
four point couplings. The principle refinement of our approach to the asymptotic dy- 
namics is that we examine the asymptotic properties of matrix elements corresponding 
to specific interactions rather than considering operators. This has the advantage of 
requiring only the machinery for the convergence of sequences of c-numbers rather than 
the more elaborate needs of operator convergence. 

In the case of 4 theory it was found that, when the incoming wave packet had test 
functions with non-overlapping supports in velocity space, the asymptotic interaction 
Hamiltonian is weakly vanishing. This condition on the test functions, of having non- 
overlapping support in velocity space, is exactly that which is required in the LSZ 
formalism and the construction of the S'-matrix. Our result is in complete agreement 
with perturbation theory and shows why it works. 

In the case of scalar QED we used our methods to show that the matrix elements 
associated with the four point interaction term are all asymptotically trivial which is 
again in line with the results of perturbation theory. For the three point interaction 
term of scalar QED, our methods show how the asymptotic dynamics associated with 
the event of a charged particle emitting a photon can be shown to be exactly the same 
as that of a charged particle with known non-trivial asymptotic dynamics, and this 
conformed to the approximation given by KF. The spin independence of this result 
immediately translates to the fermionic theory. 

What can we learn from this work about QED? Firstly that the coupling does not 
'switch off' at large times. KF further showed that this implies that the Lagrangian 
matter field does not asymptotically approach the free field of the plane wave expansion. 
Rather there is a distortion factor which expresses itself in perturbation theory in the 
branch cuts (instead of poles) in the matter field two point function. They drew 
the conclusion from this that it is not possible to describe charged particles in QED. 
Although this paper supports the non- vanishing of the interaction, we feel that this last 
conclusion is not justified. What one needs is to find the fields which do asymptotically 
approach the plane wave expansion and can therefore be interpreted as particles. (It is 
in fact clear from the start that the Lagrangian matter cannot hope to do this since it 
is not gauge invariant.) That such fields exist has been shown elsewhere ITBfl and that 



their Green's functions have a good pole structure has been amply demonstrated, see, 
e.g., [pTTj— |19H . These are the physical fields which should be identified with the charged 
particles seen in experiment. 

We would like to suggest here three further areas for study: massless QED is a 
theory with collinear divergences and as such a playground for understanding QCD. 
The asymptotic dynamics of this theory |2(J requires further study, in particular the 
physical asymptotic fields need to be constructed. Finite temperature field theory is 
another area where infra-red divergences are important, here, of course, the residual 
asymptotic dynamics of zero temperature will be acerbated by excitations from the 
heat bath. Finally in QCD confinement shows that the interaction does not switch 
off and the strong interaction between quarks and gluons is indeed supposed to grow 
with the separation. The application of the methods of asymptotic dynamics and the 
construction of physical fields at short distances fl2l] , |22l| could have implications for jets 
production. It has though been demonstrated that there is a topological obstruction to 
the construction of an isolated quark or gluon [22|, how this relates to non-perturbative 
effects in the asymptotic dynamics of QCD is a topic for future work. 
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A Appendix 

In this section we shall provide a proof of Theorem [I]. This is split into two parts: First 
we shall show that the two integrals 



d 3 qd 3 kf(q)g(q)q' 1 h^k)e- i,p ' t and j d 3 qd 3 kf(q + k)g(q)q' 1 h^k)e- i ' p ' t 
are asymptotically equivalent. Then we shall show that the two integrals 

d 3 qd 3 kf(q + k)g(q)q»h^k)e-^' t and / ' d 3 qd 3 kf(q + k)g(q)q^h il {k)e- i ^ t 



are asymptotically equivalent. Together, they give the proof of the theorem. The first 
of these results is 



Lemma 1 Let /, p,^, ^' be as defined in (pOj). Then the two integrals 

d 3 qd 3 kf(q)g(q)q tl h^k)e- i ' p ' t and J d 3 qd 3 kf{q + k)g(q)q fl h IM (k)e-^' t (32) 
have the same asymptotic limit, i.e., if It is defined by 

-iip't 



I t = J d 3 qd 3 k [f(q + k) - f(q)]g(q)q"h tl (k)e 



(33) 



then I t — > as \t\ — > oo. 

Proof We first take a fixed value of q and consider 

l t (q) = J d 3 k [f(q + k) - /(g)]g'V*)e~ Vt • ( 34 ) 

This integral is clearly well defined and there is a a positive number M, say, such that 
J d 3 k\q"h ll (k)\ < M . 

Now given e > 0, choose a 5 > such that \f(q + k) — f(q)\ < ^ if |fc| < <5 . Let 
U\ = {k : \k\ < 5} and C/ 2 = {k : |fe| > 6/2} and let pi,/?2 be a smooth partition of 
unity subordinate to U\, Ui respectively. Write 



ll(q) = J d 3 k Pl (k) (f(q + k)-f(q))q^(k)e-^ 
l 2 t (q) = J d 3 kp 2 (k) (f(q + k)-f{q))<fhr(k)e-« i 



(35) 



For the first of these we have \Xj(q)\ < e/2 by construction. For the second, the inte- 
grand in X t 2 (qr) is defined on U2, which does not contain zero, so that ip' is differentiate 
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on the domain of integration (i.e., U2) ■ One can then apply the method of stationary 
phase to this integral, using the partial derivatives of ki, to show that 2T t 2 (qr) — * as 
\t\ — > 00 . 

We have shown that T t (q) — > as \t\ — > 00, for every q. Note that T t {q) satisfies 
the inequality |X<(q)| < / d 3 k \f(q + fe) — f(q)\ |g M /i M (fc)|, and the latter is in L 1 (qr). 
The Lebesgue Dominated Convergence Theorem |23| may now be invoked to prove 
that I t —>■ as \t\ — > 00, and this completes the proof. □ 

Lemma [l] means that we have to show that the two integrals 

d 3 qd 3 kf(q + k)g(q)q' I h^k)e' l1p ' t and / d 3 q d 3 kf(q + k)g(q)q fl h^k)e-^ t (36) 



have the same asymptotic limits, where ip, ip' are defined in fl3"0|) and (|27|) respectively. 

Before proceeding with this, let us recall the form of Taylor's theorem for a smooth 
function u, i.e. 

u(q + k) - u(q) - diu(q)ki = R{q,k) 

and \R(qM^ &g |fc| ^ _ (3?) 

\k\ 

Thus, for a fixed value of q and a given £ > 0, there is a 5 > such that if \k\ < 5 
then \R(q, k)\/\k\ < e . If the choice of fc is restricted to the set {& : \k\ < 1}, then 
3i > such that V \t\ > t , we have \R(q, k/t)\/\k/t\ < e and from this we are able 
to conclude that \tR(q, k/t)\ < e for all \t\ > t . 
We shall now prove our final theorem. 

Theorem 2 Let f,g,h^,ip,ip' be defined as in and (WW- If we define I t as 



I t = J d 3 qd 3 kf(q + k)g(q)q% l (k)(e-^- e -^ t ), (38) 
then I t —>■ as \t\ — > 00. 

Proof In the following discussions the value of the functions /, g, in fl54]) are not 
important and the only property that is required of them is that they are test functions. 
For convenience, therefore, we shall write the integral as 

I t = [ d 3 q d 3 k h(q, k) (e" #t - e" # '*) . (39) 



with h(q, k) being a test function. 
Take a fixed value of q and write 



I t (q) = J d 3 k h(q, fc)e" #t (l - e^"^') 

= Jd(j \ h(q, j^j e -W*)t(i - e^ k ^) . (40) 

As in Lemma |l|, we shall first show that h{q) has a zero asymptotic limit. In the 
latter integral we have changed the variable from k to k/t so that now we are writing 
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ip{k/t) = E q+k/t - E q 



UJ 



(fe/t) 



and $(fe/t) = E q+k/t -E q -(q- k)/(tE q ) . The latter 



expression is in a form that will allow the use of Taylor's theorem. 

Let U\, U 2 be the open cover of M 3 given by U\ = {k : \k\ < 1}, U 2 = {k : \k\ > |}, 
and let pi,p 2 be a smooth partition of unity subordinate to this cover. We write the 
second integral in fl40| ) as the sum of two integrals by incorporating this partition, i.e. 



d 



Pi{k)h\q 



k 



-iip(k/t)t. 



(41) 



k 



d[-\ p 2 (k)h[q,-\(e 



k 



-iip(k/t)t 



-iip'(k/i)t 



and we shall deal with I}{q) first. 

Due to the presence of pi, the integral !}{q) is defined on {k : \k\ < 1}. Now given 



e > 0, 3 5' > such that if \6\ < 5' then |1 — e | < e . From Taylor's theorem, (see 
and the paragraph following it) we can find a to > such that, for all |fe| < 1, if \t\ > to 
then | (E q+k/t - E q - (k ■ q)/{tE q )) t\ = \tR(q,k/t)\ < 5', say. Then if \t\ > t , we 



have 



e = e 



d 3 k\h(q, k)\ , 



(42) 



where the latter integral is well defined, is independent of t and has been obtained from 
the previous one by changing the variables. 

The next step is to show that the asymptotic limit of If{q) is zero. For this we 
consider the last integral in (JT) as the difference of the two obvious integrals given by 
their exponential terms, i.e. 



k 



d \ T P2(k)h\ q, - e 



k 



-i-ip(k/t)t 



d 3 kp 2 (kt) h(q, k) e 



(43) 



I 2 2 (t,q) 



k 



d \ T P2(k)h[ q,- e 



k 



-iip'{k/t)t 



d 3 kp 2 (kt) h(q, k) e 



-iip'(k)t 



with the final forms of If, l\ being obtained by the obvious change of variables. 

We shall first examine l\ and we begin by changing the variable for k again, this 
time to k = uk with \k\ = 1 and uj > 0, i.e. polar coordinates. We now have 



A 2 (t,q) 



dk / u) 2 dio p 2 (kujt) h(q, kuot) e 



-ii/j(ku!)t 



(44) 



Note that p 2 (kt) = p 2 (kut) and this will be zero for < ut < | . The exponent of 
the integral in ([44]) is 



tjj(ku) 



E q+ku ~ E q -U 
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so that 



c , , def dip q- k + u 

Ww) = a; =_ ^r _ (45) 



As w -> 0, £i(^) tends, uniformly in k, to (q • k/E q ) — 1, and since \q ■ k/E q \ < 
(\q\/E q ) < 1, we have ^i(w) = dip/du is bounded, uniformly in k, strictly away from 
in a neighbourhood of u — . It is also easy to check that £i(u;) is non-zero for any 
finite value of to . 
We now have 

= 1 [ dk r ' diu^(p 2 Ckiut)u 2 h(q,kLu)-^-)e-^\ (46) 

with the latter expression being obtained after integration by parts in to, and noting 
that the boundary terms vanish. This can then be written as the sum of two integrals 

+ l -[dk r du, f ( ^MilM' ] e -*(M< . (47) 



it 7 Jo \daj J \ ^i(cu) 



If vi(q,ku) = u) 2 h(q, ku)/£i(u) then dv^q, koS)/du is rapidly decreasing at in- 



finity. The first integral in ( f4?[ 

(l/|t|)/dfc fo°°duj dvi(q,ku>)/du> 



can now be disposed of since it is bounded by 

and since this is well defined, vanishes as \t\ — > oo . 

Before dealing with the second integral in (f|7|) it is worthwhile examining the prop- 
erties of the derivatives of p%. We can write this function as p 2 (kuj) = p(k,u), empha- 
sising the fact that p2 is a function of two variables. Now let us write 

g{ku) = ^-p 2 (kuj) = d 2 p(k, w) (48) 

Since p 2 is smooth and p 2 (ku) = 1 for u > 1, the function g(ku) is a smooth function 
which has the important property that it is zero for oj > 1, i.e., g has its support in 
{oj : 1/2 < < 1}. 

The second integral in (f4~T|), which we denote by J~i(q), is now 

= -J dk J uj 2 duj—(p 2 Ckujt)^ Vl (q,kuj)e- l ^ w)t 

= L f die f u 2 dcjtofat)v 1 (q,iw)e* l '$ uyt (49) 



o 



OC 



- I dk I u) 2 du gCktot) Vl (q,kuj)e' ilp{kuj)t 



o 



The last integral in fl4~9"D is now written in a form which is more convenient and which is 
obtained by the following changes of variables. First replace kuj by k, and then replace 
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k by k/t . The form of J{{q) then changes to 

^) = -J d (j) 7) (50) 

In order to deal with (|50|) we shall first have to examine the second integral in (^) . 
Applying the same methods to J| as we have to If one can easily show that if can also 
be written as the sum of two integrals, as in (f47|), but with the exponent if) replaced 
with if)' and £1 replaced with £ 2 > where 

c f v def 0V' ,f \ q ~ k 

= = -g — ( 51 ) 

Thus 

+i I A jTu^|-(,, 2 (M) fe^) e-*«"", (52) 

c.f. (^7|). The function £ 2 has similar properties to £1 and the first integral in ( |52"D 
vanishes asymptotically in a similar fashion to the corresponding integral in (f47|). This 
leaves us with the latter integral in (|52"D, which we shall write as J^{q) ■ Using the 
same changes of variables, this can can be written as 

J*(q) = IJ<? e (k) v 2 (q, e-^» (53) 

The final objective is to show that Jl(q) — <J£(q) is asymptotically vanishing and 
this can be achieved in two steps. 
First define J£{q) as 

JHq) = l -jd" g(k) Vl (q, ^ e'^ . (54) 

i.e., the v% in J7 2 *(q) is replaced by V\. Then it is straightforward to prove, along the 
lines of Lemma [L], that, for all q, we have J\{q) — <J£(q) — ► as \t\ — > 00. 

Secondly, we must show that Jl(q) — >J£(q) ^ as \t\ — ► 00. This is also straight- 
forward and can be fashioned along the lines of the proof that I}(q) — > as \t\ — > 00 
(see (PI) and the paragraph that follows it.) We omit the details. 

This proves that, for every q, It(q) — > as \t\ — > 00 (see(f4"0"|)). The Lebesgue 
Dominated Convergence Theorem now implies that It — > as \t\ — > 00, as required. □ 
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